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Abstract. The main objective of this article is to study the effect of the 
moisture on the planetary scale atmospheric circulation over the tropics. The 
modeling we adopt is the Boussinesq equations coupled with a diffusive equa- 
tion of humidity and the humidity dependent heat source is modeled by a 
linear approximation of the humidity. The rigorous mathematical analysis is 
carried out using the dynamic transition theory. In particular, we obtain the 
same types of transitions and hence the scenario of the El Nifio mechanism 
as described in |10l . The effect of the moisture only lowers slightly the 
magnitude of the critical thermal Rayleigh number. 
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1. Introduction 

This article is part of a research program to study low frequency variability of 
the atmospheric and oceanic flows. As we know, typical sources of climate low fre- 
quency variability include the wind-driven (horizontal) and thermohaline (vertical) 
circulations (THC) of the ocean, and the El Nirio Southern Oscillation (ENSO). 
Their variability, independently and interactively, may play a significant role in 
climate change, past and future. The primary goal of our study is to document, 
through careful theoretical and numerical studies, the presence of climate low fre- 
quency variability, to verify the robustness of this variability's characteristics to 
changes in model parameters, and to help explain its physical mechanisms. The 
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thorough understanding of the variability is a chaUenging problem with important 
practical implications for geophysical efforts to quantify predictability, analyze error 
growth in dynamical models, and develop efficient forecast methods. 

ENSO is one of the strongest interannual climate variability associated with 
strong atmosphere-ocean coupling, with significant impacts on global climate. ENSO 
consists of warm events (El Nino phase) and cold events (La Niiia phase) as 
observed by the equatorial eastern Pacific SST anomalies, which are associated 
with persistent weakening or strengthening in the trade winds; see among others 
[IllllillllllllinilllllllllllHlIinillSllSlin]- An interesting current debate 
is whether ENSO is best modeled as a stochastic or chaotic system - linear and 
noise-forced, or nonlinear oscillatory and unstable system [E]? It is obvious that a 
careful fundamental level examination of the problem is crucial. For this purpose, 
Ma and Wang JlO\ initiated a study of ENSO from the dynamical transition point of 
view, and derived in particular a new oscillation mechanism of the ENSO. Namely, 
ENSO is a self-organizing and self-excitation system, with two highly coupled os- 
cillation processes-the oscillation between metastable El Nino and La Nina and 
normal states, and the spatiotemporal oscillation of the sea-surface temperature 
(SST). 

The main objective of this article is to address the moisture effect on the low 
frequency variability associated with the ENSO. First, as our main purpose is to 
capture the patterns and general features of the large scale atmospheric circulation 
over the tropics, it is appropriate to use the Boussinesq equations coupled with 
a diffusive equation of humidity. In addition, the humidity effect is also taking 
into consideration by treating the heating source as a linear approximation of the 
humidity function. 

Second, although the introduction of the humidity effect leads to substantial dif- 
ficulty from the mathematical point of view, we haver shown in Theorem 4.1 that 
the humidity does not affect the type of dynamic transition the system undergoes. 
Namely, we show that under the idealized boundary conditions, only the continu- 
ous transition (Type-I) occurs. However, the critical thermal Rayleigh number is 
slightly smaller than that in the case without moisture factor. To see this effect of 
the humidity, we refer to formula p.3ip which reads 



Le I 24 / Le aqa\^KT 



Since in the expression of thermal critical Rayleigh number R^, the coefficient of 
the humidity Rayleigh number R is negative, and this indicates that the presence 
of humidity lowers down the critical temperature difference for the onset of the 
dynamic transition. 

Third, we remark that the perturbation analysis in _10 can be applied to the 
case here to carry out the analysis, and we can show that that under the natu- 
ral boundary condition, the underlying system with humidity effect will undergo 
a mixed type transition. In addition, as we argued in [10], it is necessary to in- 
clude turbulent friction terms in the model to obtain correct convections scales for 
the large scale tropical atmospheric circulations, leading in particular to the right 
critical temperature gradient and the length scale for the Walker circulation. 
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Finally, based on these theoretical results, it is easy then to conclude the same 
mechanism for ENSO as proposed in Ma and Wang [?]. Namely, ENSO is a self- 
organizing and self-excitation system, with two highly coupled oscillatory processes: 
1) the oscillation between the two metastable warm (El Nino phase) and cold events 
(La Nina phase), and 2) the spatiotemporal oscillation of the sea surface tempera- 
ture (SST) field. The interplay between these two processes gives rises the climate 
variability associated with the ENSO, leads to both the random and determinis- 
tic features of the ENSO, and defines a new natural feedback mechanism, which 
drives the sporadic oscillation of the ENSO. With the deterministic model consid- 
ered in this article, the randomness is closely related to the uncertainty/fluctuations 
of the initial data between the narrow basins of attractions of the corresponding 
metastable events, and the deterministic feature is represented by a deterministic 
coupled atmospheric and oceanic model predicting the basins of attraction and the 
sea-surface temperature (SST). 

It is worth mentioning that as practiced by many in the field, the randomness 
associated with the ENSO may be incorporated by introducing the noise directly 
in the model. In addition, from the predictability and prediction point of view, it 
crucial to capture more detailed information on the delay feedback mechanism of 
the SST. For this purpose, the study of explicit multi-scale coupling mechanism to 
the ocean is inevitable. 

This article is organized as follows. Section 2 gives the objective Boussinesq 
model with humidity. The eigenvalue problem is analyzed in Section 3. The transi- 
tion theorem is stated and proved in Section 4. In section 5, the turbulence friction 
factors are considered and the corresponding critical temperature difference and 
the wave numbers of Walker circulation are checked. 

2. Model for Atmospheric Motion with Humidity 

2.1. Atmospheric circulation model. The hydrodynamical equations governing 
the atmospheric circulation is the Navier-Stokes equations with the Coriolis force 
generated by the earth's rotation, coupled with the first law of thermodynamics. 

Let {ip, 6, r) be the spheric coordinates, where ip represents the longitude, 9 the 
latitude, and r the radial coordinate. The unknown functions include the velocity 
field u = {ui^,ue,Ur), the temperature function T, the humidity function q, the 
pressure p and the density function p. Then the equations governing the motion 
and states of the atmosphere consist of the momentum equation, the continuity 
equation, the first law of thermodynamic, the diffusion equation for humidity, and 
the equation of state (for ideal gas) which read 
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Here < f < 2tt, —tt/2 < 9 < tt/2, a < r < a + h,a is the radius of the earth, h 
is the height of the troposphere, is the earth's rotating angular velocity, g is the 
gravitative constant, /i, kt, Kg, Cy, R are constants, Q and S are heat and humidity 
sources, and k — (0, 0, 1). The differential operators used are as follows: 

1. The gradient and divergence operators are given by: 

f d ld_ d_ 

r cos 9 d(p ^ r d9^ dr 

1 9 , o , 1 d(u0cos9) 1 duu, 

div u = -^-^{r^Ur) + TTT. + H o ; 

or r cos 9 09 r cos 9 dip 

2. In the spherical geometry, although the Laplacian for a scalar is different 
from the Laplacian for a vectorial function, we use the same notation A for both 
of them: 

/ 2 dur 2 sin 9 dug Uu, 

Au ^ [Au^ + — -— + — — — , 

V cos 9 0(p r^ cos^ 9 dip cos^ 9 



Aug 
AUr 



2 dur ug 2 sin 9 du 



j,2 QQ J.2 (2os2 cos^ 9 dip ' 

2ur 2 d{ug COS 9) 2 du^ \ 

COS 9 d9 cos 9 dip / ' 

1 df , ^d , 1 d^f 1 9/, 2 5 ^ 

cos 9 d9 d9 cos^ 9 dip'' r^ dr dr 

3. The convection terms are given by 

/ U^Ur U^Ug 

VuU ^ [u ■ vuip H — tan 9, 

\ r r 



UgUr ^1 ul+l 

u ■ wug H 1 — — tan 9, u ■ VUj, 

r r r 

4. The Coriolis term 2f1 x u is given by 

2fl X u — 2fl(cos 9ur — sm9ug, sin^u^, — cos^Ui^). 

Here (l is the angular velocity vector of the earth, and ft is the magnitude of the 
angular velocity. 

The above system of equations were basically the equations used by L. F. Richard- 
son in his pioneering work |16j . However, they are in general too complicated to 
conduct theoretical analysis. As practiced by the earlier researchers such as J. Char- 
ney, and from the lessons learned by the failure of Richardson's pioneering work, 
one tries to be satisfied with simplified models approximating the actual motions 
to a greater or less degree instead of attempting to deal with the atmosphere in all 
its complexity. By starting with models incorporating only what are thought to be 
the most important of atmospheric infiuences, and by gradually bringing in oth- 
ers, one is able to proceed inductively and thereby to avoid the pitfalls inevitably 
encountered when a great many poorly understood factors are introduced all at 
once. The simplifications are usually done by taking into consideration of some 
main characterizations of the large-scale atmosphere. One such characterization is 
the small aspect ratio between the vertical and horizontal scales, leading to hydro- 
static equation replacing the vertical momentum equation. The resulting system 
of equation are called the primitive equations (PEs); see among others [7]. The 



TROPICAL ATMOSPHERIC CIRCULATIONS WITH HUMIDITY EFFECTS 



5 



another characterization of the large scale motion is the fast rotation of the earth, 
leading to the celebrated quasi-geostrophic equations [2]- 

2.2. Tropical atmospheric circulation model. In this article, our main focus 
is on the formation and transitions of the general circulation patterns. For this 
purpose, the approximations we adopt involves the following components. 

First, we often use Boussinesq assumption, where the density is treated as a 
constant except in the buoyancy term and in the equation of state. 

Second, because the air is generally not incompressible, we do not use the equa- 
tion of state for ideal gas, rather, we use the following empirical formula, which can 
be regarded as the linear approximation of (12.51) : 

(2.6) p = po[l - ariT - To) - ag{q - qo)], 

where po is the density at T = Tq and q — qo, and ar and Uq are the coefficients 
of thermal and humidity expansion. 

Third, since the aspect ratio between the vertical scale and the horizontal scale is 
small, the spheric sheU the air occupies is treated as a product space S"^ x (a, a + h) 
with the product metric 

ds^ = a^de^ + sin^ Bd(^^ + dz^ 

This approximation is extensively adopted in geophysical fluid dynamics. 

Fourth, the hydrodynamic equations governing the atmospheric circulation over 
the tropical zone are the Navier- Stokes equations coupled with the first law of 
thermodynamics and the diffusion equation of the humidity. These equations are 
restricted on the lower latitude region where the meridianal velocity component uq 
is zero. 

Let (0, z) G M — (0, 27r) x (a, a + /i) be the coordinate where is the longitude, a 
the radius of the earth, and h the height of the troposphere. The unknown functions 
include the velocity field u = {u^,Uz), the temperature function T, the humidity 
function q and the pressure p. Then the equations governing the motion and states 
of the atmosphere read 



9u0 
dt 
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dt 
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2 du^ 

a 0(p 
- ariT - To) - aq{q - qo)), 


2m2, 


- 2Qu^ — 



1 dp 

poa dcj)' 
]_dp 
Po dz 



(2-7) dT 

— + {u-V)T^ktAT + Q, 
dt 

-^ + {u-\7)q^ KqAq + H, 

1 du^ ^ duz_ _ Q 
a d4> dz ' 

where Q, is the earth's rotating angular velocity, g the gravitational constant, po the 
density of air at T = Tq , v the viscosity, a-r and aq are the coefficients of thermal 
and humidity expansion, kt and Kq the thermal and humidity diffusion coefficients, 
and Q, H the heat and humidity sources. However, the differential operators used 
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in (|2.7|) are as follows 



"0 
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d 
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2 du:, 2 2u^ 

where A/ is the Laplacian for scalar functions, and Au is the Laplace-Beltrami 
operator on (0, 27r) x (a, a + h) with product metric. Based on the thermodynamics, 
the heat source Q is a function of the humidity q. For simplicity, we take the linear 
approximation 

(2.8) Q^ao + aiq, 
and the humidity source is taken as zero 

(2.9) H ^0. 

The boundary conditions are periodic in (/)-direction 

(2.10) [u, T, q){(t> + 27r, z) = (w, T, q){cp, z), 

and free-slip on the earth surface z = a and the tropopause z — a -\- 

Uz = 0, =0, T ^Tq, q = qo, at z = a, 

(2.11) 

Uz=0, ^=0, T = Ti, g = 0, atz = a + /i, 
oz 

where Tq, Ti, qo and qi are constants satisfying 

(2.12) To>Ti, qo>0. 

The problem (|2.7p - (|2.1ip possesses a steady state solution 
u = 0, 

^ 7l 3 7o 2 , 

T = - — z-" - - — z^ + ciz + Co, 

(2.13) ~ ?o , go., , X 

P^^J Pog{^-aT{T -To)-aq[q-qa)^dz, 

where 

70 ao + ai(7o(^ + a)//i, 7i = ai^o/^, 

Co = - - — 0.-' + - — a +Ta\--[-- — [a + hY + - — [a + h) + Ti 



ci = -Uto - Ti) - — (^{h^ + 3ha + 3a^) - ^{h + 2a) 
n kt ^ o 2 
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To obtain the nondimensional form, let 

X = hx' , a = hro, 

t = h^t' /kt, u = ktu' /h, 

T={To-Ti)T' + f, q = qoq' + q, 
p = povktp' /h^ +p, 

and consider 

{xi,X2) = iro(f>,z), (ui,U2) = [u^,u^). 
Omitting the primes, the equations (|2.7p become 



dui ( 2 du2 1 dp\ 
— - = Pr Aui H - Tr~ 

Ot V Tq OXi Tq OXi / 

— ijJU2 — [u ■ V)wi U1U2, 

ro 

du2 p / . 2 dui 2 . „^ . p dp\ 

— - = Pr Au2 2 "2 + RT + Rq- -— 

ot V ro 0x1 Tq 0x2 

(2.14) +a;ui - (u • V)w2 + — u?, 

ro 

^ = AT - U2 + aq-{u- V)T, 

at iq — li ax2 

dq 

— =Le Aq + U2 - {u- V)(j, 

dui _^ du2 _ Q 
dx\ dx2 

where the nondimensional physical parameters are 

V 

Pr = — the Prandtl number, 

R = " ^ — the thermal Rayleigh number, 

R = — the humidity Rayleigh number, 



(2.15) 



Is. 



Le — — - the Lewis number, 



2nh'^ 

Lu — the earth rotation, 

KT 

aiqoh? 



kt{To - Ti) 



and 



(2.16) 



1 df{hx2) 7ofe' , 1, 



Tq-Ti dx2 kt{To-Ti) 2 

amh"^ ,22 In 



2KT{To-Tiy ^ " " 3^ 

where ro < X2 < rg + 1, and ro is the nondimensional radius of the earth. For 
simplicity, we take the average approximation X2 — + ^. In this case, (|2.16p is 
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given by 



1 



dT 



1 



To - Ti dX2 X2=ro + ^ 

Thus, the equation (j2.14p is approximated by 

dui 
'dt 



(2.17) 



Pr Aui H :rui - - — 

rg OXl Tg OXi J 



■ UJU2 — {u ■ V)ui 1*1^2, 
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dt 



dT 
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=Prf Au2 - — ^ - + RT + Rq 

V ro dxi 

1 2 

+ UJUl - (U ■ V)U2 H Ui , 

=AT + 71*2 + aq-{u- V)r, 
=Lc Aq + U2- {u- V)(j, 



dt 

div u — 



dx2' 



where 
(2.18) 



7 = IH . 

^ 24 



The domain is M = [0, 27rro] x (ro, j'o + 1), and the boundary conditions are given 

by 

{u,T,q){xi + 27rro,X2) = (u, T, g)(a;i, 2:2), 

(2.19) dui 

U2 = 0, - — = 0, T = 0, g = 0, at X2 = ro, ro + 1. 

dX2 



For the problem p.l7p - p.l9p . we set the spaces 

H = |(u,r,g) e L^{M,M.'^) div u = 0, U2 = at ro,ro + l|. 



(2.20) 



Hi = \ {u,T,q) e H^{M,R^)nH {u,T,q) satisfies (1^1^ | 



3. Eigenvalue Problem and Principle of Exchange of Stability 



3.1. Eigenvalue Problem. To study the transition of (|2.17p - (|2.19p from the basic 
state, we need to consider the fohowing eigenvalue problem 



(3.1) 



„ / . 2 du2 

V ro 0x1 



2 9ui 2 „^ n 9p \ 

' fr n 



Prf Au2 TT"^ 9'"2 + RT + Rq — -^-^ ) + ojui — /3u2, 

V ro dxi Tq dx2 ' 



AT + 71*2 aq — /3T, 
Le Ag + U2 = ;Sg, 
9ui 9u2 



i9a:i dx2 



= 0, 
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supplemented with the boundary conditions ()2.19p . By the boundary conditions 
()2.17|) . the eigenvalues of can be solved by separation of variables as follows 



(3.2) 



(3.3) 



^2 
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, N . kX\ 

= 9fc(a;2)sin , 

ro 





p =^fc(x2)sin \-Bk{x2)cos 

ro ro 

Based on the continuity equation in (13. ip . we have 

ro dvk 



(3.4) 
Let 
(3.5) 



k 'dx2 



ujro 2 

= TTr'"k{X2), Ak=Pk^ Vk{X2)- 

KFr ro 



Due to (I2.19p . plugging p.2p or p.3p into p.ip . we obtain the following system of 
ordinary differential equations 

1 



(3.6) 



Pr D^Uk K'^k Pk = (iuk, 

V r^ ro / 

¥r(Dlvk - \vk + RTk + Rqk - DpA ^ /Svk 



DlTk + "fVk + aqk = /3Tk, 
he Dlqk +Vk ^ Pqk, 

Duk =0, Wfc = 0, Tk =0, qk^ 0, at X2 = ro, ro + 1, 



where 
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dX2 
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Plugging Vk = sin j7r(x2 — tq) into p.6p . we see that the eigenvalue /3 satisfies the 
cubic equation 

("fc, + /5)(Lc al^ + P){PTal^ + ^ + P)al^ + ^^{a% + /3)(Le a% + /3) 
(3 7) ° ° 

^(7Le a^^. + 7/3 + a) —{a% + /?) = 

where 

(3.8) al^. =/vr2 + ^ (j>l,fc>l). 

3.2. Principle of exchange of stability. The linear stability of the problem 
(I2.17l) - (|2.19p is dictated precisely by the eigenvalues of p.ip which are determined 
by p.7p . The expansion of p.7p is 

/33 + [(1 + Le + Pr)a2 +^ + 

Pr Prfc^ 
+ [(Pr + Le Pr + Le)^^,. + -j(l + Le)(<.2^. + -2 ) - -j;^{lR + i?)]/? 
(3.9j '0 '0 'o"fei 

Pr Prfc^ ~ Le at 
+ Le atj(Pra?,^ + — ) + -^^{~a%R -aR^ 7Le a^^i? + ^) 

= 0. 

We see that /3 = is a solution of p.9p if and only if 

Pr Prfc'^ Le at 
(3.10) Le afej(Pra^.^ + — ) + -^r^{-a%R - aR - jLe al^R + j^) = 0. 

Inferring from (I^TTS)) and (PT5)) . 

aiargqoh^ aiarh'^ 
a/t = ^ = R, 

7i? = i?+"™^' Wfe ^ 

Hence, we rewrite p.lOp as 

We define the critical Rayleigh number Rc as 
(3.13) 

7-: + (l+ )7T 2 M? + «fe7T^(" 



i?c = min 

fej>i 



2^ ' „2 



: mm 

;d>i 



Vl^ ^ ^ "^^Le a,a%nT)'' "fe. p ^"fe. - " ,^ 

(t^ + (1 + ^) "^""2"' )^ + «ti|(«L + ^) + 4 
VLe 24 ^Leaga^^KT/ rl' rl 



where 

_2 



'0 

Based on - ([XT^ . by Definition ([XTO)) . we shall prove the following PES later. 
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Lemma 3.1. Let I3jk he the eigenvalues of LS. 1]) that satisfy L3.9\) . Let kc > 1 be 
integer minimizing 113.13]) . and I3\ (1 < i < 3) be solution of 113. 9\) with {k,j) — 
{kc, I), and 

mc(3l<mtpl<D\zi3l. 

Then /3^. must be real near R — Rc, and 



< 0, R< Rc 
= 0, R^ Rc 
> 0, R> Rc 



Moreover, we have 



UicPi (Rc) <0, for J ^2,3. 



Pjk{Rc) < for Pjk being real. 



Remark 3.1. The value Rc defined by (|3.13p is called the critical Rayleigh number 
of real PES, which provides the critical temperature difference ATc = To — Ti by 
R = Rc, i.e., 

axgh^l^Tc 



Equivalently, we have 



KTl' 



(3.14) 



Rc- 



24 



Le aqU'l^^KT'' 



,4 '0 (2 , \, "fed 



where fee > 1 is the integer which satisfies p.l3p . 

Next, we consider the PES for the complex eigenvalues of (|3.ip . Let (3 = ipo 
{fa 7^ 0) be a zero of ([231). Then 

Pr k^ Prfc^ 

pI =(Pr + Lc Pr + Lc)al^ + —(1 + Le){al^ + — ) " ZJ::^^^^ + -^)' 



pI 



Pr, 



Prfc2 



Le at 



Pr 

(1 + Le + Pr)Q;^j + — 



Hence, Equation p.9p has a pair of purely imaginary solution ±ipo if and only if 
the following condition holds true 

Pri 



(1 + Lc + Pr)afej- 



(3.15) 



+ 



Pr 



(l + Le)(aL + -2) 



'0 
2Pr. 



(Pr + Lc + Le VT)alj 
Prfc2 ^ 

Prfc2 



Le aljiPialj + — ) + -^-^{-aljR - aR - -/Le af^R - 



'0 



> 0. 



'0 
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It follows from ([3lT|) and (|3?T5|) that 



(3.16) 



Pr 

— + (Pr + l)a: 



R 



aiaxh? ( 1 /Pr ,^ , , 
QfgKT V 24Vro 



/Pr 



(Pr + Le)a2, 



1|(M1 + Lc, 



(1 + Lc) 



Pr 
Pr/c^ 



Pri 



(1 + Lc + Pr)a^^. + — (1 + Le)a^^. 



s 9 Pr 



where af^ is as defined in (|37 

According to p.l6p , we define the critical Rayleigh number for the complex PES 
as follows 



K 



mm 

fc,j>i 



^ + (Pr+l)a2 



■X 



(3.17) 



r \a\ctTh'' ( i /Fr 



kj 

1 /Pr 



l)a^ 



+ 



+ 



OiqKT 
1 9 

fc2pr 



(1 + Le) 



+ Le) 







— + (Pr + Le)a2, 



i? 



fc^Le 



Pr 
Pr/fc2 



(1 + Lc + Pr)«2^. 



Pr 
72" 



(1 + Le)at, 



(1 



Pr 
'0 



)}. 



Then we have the following complex PES. 



Lemma 3.2. Let k*, j* be the integers satisfying ^3.17^ , and and Pf,,-, be 

the pair of complex eigenvalues of \3. 9\l with {k,j) = ik*,j*) near R — R*. Then, 



<0, R<R*, 
— 0, R — i?*, 
>0, R>R*.. 



Furthermore, for all complex eigenvalues Pkj of (13. ip . we have 
5ne/3fc,(i?:)<0, for(fc,j)^(fc:,j:). 



3.3. Proof of Lemmas 13.11 and 13. 2i We notice that all eigenvalues of p.ip are 
determined by (j3.9p and the eigenvalue equations 

AT + = /3r, 
(3-18) , / „ 

^ ' Le Ag = (3q, 

with the boundary condition (|2.19p . It is clear that the eigenvalues of p.lSp are 
real and negative. Hence, it suffices to prove Lemmas 13.11 and 13.21 for eigenvalues 
hj of (131]). We rewrite (|3H) as 

(3.19) ^3 + 62^2 + 61/3 + 60 = 0, 
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where 

Pr fc^Pr 

'o 'o"fej 

Pr fc^ PrA:^ 

bi =(Pr + Le Pr + Le)a|j + —{1 + Le){alj + — ) - -^-^{-/R + R), 

Pr Prk^ Le ctf ■ 
bo =Le atjiPTalj + — ) + ^-^(-"fej-R - ai? - 7Le a^^-i? + 2"^). 

We see that 

(3.20) bo, 6i, &2 > at i? = for all fc, j > 1, 

which implies that all real eigenvalues of p.9p are negative at i? = 0. By the 
definition of Rc, we find 

6o>0,V(fc,i)7^(fc„l), 0<R<Ra, 

> 0, R< Rc 

^^•^^^ bo I =0, = for (fc,j)-(fc„l). 

< 0, i? > Rc, 

Since 6o = -PljPljPlp Lemma O follows from (15:^ and (IX^ . 
Next, let Pk'j*, the solution of p.9p . take the form 

/3,;,;(i?) = A(i?)+Zp(i?), 

^' A(i?) ^ 0, p(i?) ^ po, asR^R*. 

Inserting ((3?22)) into ((3?T9| . we obtain 

(-3p2 + 6i)A + bo- &2P' + o(A) = 0, 
-p^ + pbi + 2pb2\ + o(A) = 0. 
Since po 7^ 0, we infer from p.23p that 



(3.23) 3 



wpx^ ^0 - bo-hb2-2qX 

MR) + 0(A) = ^^^^ = + 6b2A + 

Thus, we have 

(3.24) A(i^)= ^°:,^^^^^-f' +o(A). 

2oi + 602A 

Under the condition (|3?T5l) . we see 61 > 0. It follows from (|X22l) and (|3?24| that 

{< 0, 5o < 6162, 
= 0, feo = ^1^2, 
> 0, bo> bib2, 

for i? near i?* . By the definition of i?* , it is easy to see that 

r <b,b2, R<R*c, 

(3.26) 5o< -6162, = 

This proves Lemma 
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By Lemmas 13.11 and 13. 2| we immediately obtain the following theorem which 
provides a criterion to determine the equilibrium and the spatiotemporal oscillation 
transitions. 

Theorem 3.1. Let Rc and R* be the parameters defined by i3.13]) and Jg. J7p 

respectively. Then the following assertions hold true. 

i) When R^ < R*, the first critical-crossing eigenvalue of the problem 113.1]) is 
P]^ given by Lemma \S.l{ i.e., 

r <o, 

(3.27) Pl{R)\ =0, 

[ >0, 

(3.28) lKe/3(i?c) < 0, 

for all other eigenvalues l3 of i3.1\) . 

ii) When R* < Rc, the first critical- crossing eigenvalues are the pair of complex 
eigenvalues P].,^-, and Pf.,^, given by Lemma \3.2\ namely, 

<0, ifR<R*, 

(3.29) mePl^,{R) = mzPl^,{R)' 



if R< Rc, 
ifR = Rc, 
if R> Rc, 



= 0, 
> 0, 



^fR^K, 
ifR>R:, 



(3.30) 5^e^(^c) < 0, 

for all other eigenvalues /3{R) of 113.1 



Remark 3.2. In the atmospheric science, the integer j* in p.29p is 1. Hence, the 
critical Rayleigh numbers Rc and R* are given by 

,Le 24 Le aqa^ ht ^ 



(3.31) Rc 

(3.32) R*c 



R 



A i2 



.("fee 



1 ^ a 



1 



^ + (Pr+l)a2. 



{ (l " ^ + (Pr + Da^O ) - + + L^)"- ) 

LL a„KT V 24\rn " J J \ri <=/ 



1 

24 



'Pr 



Pr 



R 



Pr 



riat, / 

(Le(l + Le)a^. + Pr (1 + Lc + Pr)a^^ + ^ (1 + Le^- 
(1 + Le) , 



Prfc! 



(1 + Pr)a2 



Pr 
75" 



where 
(3.33) 



— S- and 



4. Transition Theorem 

Inferring from Theorem 13.11 the system ()2.17p - p.l9p have a transition to equi- 
libria at R = Rc provided Rc < R*, and have a transition to spatiotemporal 
oscillation at R — R^ provided R* < Rc. 

Theorem 4.1. For the problem ^2.17\ l- (2.19\) we have the following assertions. 
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(1) When R < imii{Rc, R*}, the equilibrium solution (u,T,q) — is stable in 
H. 

(2) If Rc < R*c, then this problem has a continuous transition at R = Rc, and 
bifurcates from {{u,T,g),R) = {0,Rc) to an attractor = on R > Rc 
which is a cycle of steady state solutions. 

(3) As R* < Rc, the problem has a transition at R = R*, which is either of 
continuous type or of jump type, and it transits to a spatiotemporal oscil- 
lation solution. In particular, if the transition is continuous, then there 
is an attractor of i- dimensional homological sphere is bifurcated from 
((m, T, g), i?) — (0, i?*) on R > i?*, which contains no steady state solu- 
tions. 

Proof We shall prove this theorem by several steps. 

Step 1. Let H and Hi be the spaces defined by ()2.20p . We define the operators 
Lr = A + Br: Hi ^ H axiAG : Hi ^ H hy 

A^ = P( Pr( Ami + - , M^u, -l^),AT,Le Aq) , 

/ Pr 2 

(4.1) i?_R^ -^i 2""i ^ ^^"2, -Pr(— U2 - PT - Rq) + cJMi, 7^2 + aq, U2 

G(V') = -P({u ■ V)ui + —uiU2, {u ■ V)u2 - —uj, {u ■ V)T, {u ■ V)q) 
^ rQ ro 

where P : L^(M, M^) — > _ff is the Leray projection and 

4,^{u,T,q)eHi. 
Thus, the problem (|2.17l) - (|2.19p is expressed in form of 

(4.2) ^ = Lfl^ + G(V'), 

and the eigenvalue problem (13. ip with the condition (|2.19p is rewritten as 

(4.3) LflV = 13^- 



Step 2. We shall calculate the center manifold reduction for (14.21) in this step. Let 
and ?/^^ be the eigenfunctions of (14. 3p corresponding to [R), where is 
the eigenvalue of (14. 3p in the case of p.27p . Denote the conjugate eigenfunctions 
of and V'fc^ by i)],* and il^l* , i.e., 

(4.4) L*j,^i; = pui; , 1 = 1,2. 

The corresponding equations of (14.41) read 

V To oxi ri 0x2 ' 



2! 



(4.5) 



AT* +Pri?u^ = /3T*, 

Le Aq* + Pr^u^ + aT* = Pq* , 

dul ^ ^^2 _ Q 
dxi dx2 ' 
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Express € Hi as 
(4.6) 



ij^xtljl+yijl+^{x,y,R), 



where {x,y) E and ^{x,y,R) is the center manifold function of (|4.2p near 
R — Rc- Then, the reduction equations of (|4.2p are 



(4.7) 



dx 
dy 



1 



-(G(^),V'l*), 



Step 3: Computation of tpl. and * (i = 1,2) Based on p.2p - ()3.8p . we can derive 



-Uk^ cos7r(a;2 - tq) sm ■ 



(4.8) 



(4.9) 











= Vk^ sin7r(x2 - 


k^Xi 
'ro)cos , 

'"0 


n 


= Tfe^ sin7r(a;2 - 


- To ) cos 

ro 




= qk^ sin7r(x2 - 


ro) cos , 

^■o 


ul 


= cos t:{x2 


, J\i(^X\ 

— ro) cos 

ro 


ul 


= Vk^ sin7r(a;2 


, , kf^x^ 

- ro j sm 

ro 




= sin7r(a;2 


. kcXi 
- ro) sm 




= qk^ sin7r(a;2 - 


X . kcX\ 

- ro) sm , 

ro 



where 



(4.10) 



rfc„ = 



a + 7(Le < + J 



Similarly, we can also derive from (I4.5P the conjugate eigenfunctions as follows 



1, 



1 



(4.11) 



-u^^ C0S7r(x2 — ro) sin 



U2 = v'k sin7r(a;2 — ro) cos 
Ti* — T^^ sin7r(x2 — ro) cos 
q^* = ql^ sin7r(x2 — ro) cos 



ro 
kcX\ 

ro 

kcXi 

ro 
kcXi 

ro 
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(4.12) 



ul 


= u^^ COS i:{x2 


-ro] 


cos 

''0 


2* 


= v^. sin7r(a;2 


-ro) 


sm , 


212* 


= sin7r(a;2 


-ra 


kcXi 

sm 

ro 


q 


= ql^ sin7r(a;2 - 


-ro) 


kcXi 

sm , 

ro 



where 



(4.13) 



kr 



T,! = 



Pri? 



1 ' 



aPri? - Pri?(a^_^ + J 



Computation of the center manifold Junction. The nonhnear operator G 
defined in (14.11) is bilinear, which can be expressed as 



(4.14) 



dT 
dq 



where ip = (ui,U2,T, q) and i 
(4.15) 

,ro7r^ . 2kxi tt . 



p 



(ui , W2, T, <?)■ Direct calculation shows 

2 



2fc 



ro 2 



, — sin27r(x2 - ro) 



roTT 
4fc2 



[l + cos27r(a;2 -ro)],0,0) 



+ ( - — sin27r(a;2 - ro) sm — — , cos27r(a;2 - ro) cos , 0,0j 



4fc 



ro 4fc2 



ro 



(0 



4fc2 



cos ■ 



ro 



2 sin27r(a;2 - ro), ^ sin27r(a;2 - ro)) 



Since the first two terms on the right hand side of ()4.15p are gradient fields, we 

obtain 

(4.16) 

^2 2kxi TfeTT . 



G(^i,^i.) = (0,- 



' 4fc2 



ro 



2 ■sin27r(x2 -ro),-^sin27r(j;2 -ro))* 



Similarly, we can derive 
(4.17) 

1 2n / roTT^ X roTT^ 2kxi st 

G(V'fe,'0fc) = ( - -^cos27r(a;2 - ro) + — sin27r(a;2 - ro), sm — — , 0, 0) , 



4fc2 

(-^cos27r(a;2 - ro) - — sin27r(x2 - ro), — j-^sm— — ,0, Oj , 



4fc 



4fc2 



ro 



roTT^ 2kxi TfcTT . ^ 9fe7r . xt 
cos , — sin27r(a;2 - ro), — sin27r(a;2 -ro)j , 



4fc2 



ro 
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and 
(4.18) 



roTT^ 2kxi Ti, Ti qk*n 



Gitpl,ipl ) = (0, -^^cos— — , ^sin27r(x2 - To), ^ sin27r(a;2 - ro)) , 

G{ipk,'^k ) = { - ^^cos27r(a;2 - ^o) + — siii27r(a;2 - fn), -^j^" ^— , 0, 0) , 

(^Wk^Wk ) = (^^cos27r(x2 -ro) - — sin27r(x2 - ro), sm — — , 0, 0) , 

r^r i2 ,2*N /'n ''0^^ Tk*Tr qk*7T .t 

(^Wk^Vk ) = (0' TTT ' ^sin27r(a;2 -To), — sin27r(x2 - roj) , 

^ 4k^ ro 2 2 ' 



By (HT^ - (|iTfl) . we have 

(4.19) <G(V'^SVfc^),Vr* >=0, 

for 11,12,13 = 1,2. Moreover, direct calculation shows 

< G(0i,02),03 >= - < G(0i,03),02 >, 

(4.20) 

<G(0i,^^),^^ >=0, 



for (j)2, 4>3 £ i^i and i = 1,2. Since the center manifold function $(x, y) = 
0(|x|2 + |y|2), by (I41B - (I4:201) . we obtain 

-.T < G«,VL*),<i> > -y < G{^l,^i;),^> 
+ y<G{<^,^l),i^i;>+o{x^+y^), 
-x< G«,Vi*),$ > -y < G(Vi,V'i*),$ > 
+ X < G($, VfeJ, V-L* > +o(a:' + ?/')• 



(4.21) 



<G(V'),^L* >= 



Let the center manifold function be denoted by 



(4.22) '^^ '^Pli^^y)'l'lr 



By (Oa - g^H), only 

^0^2 =(0,0,0,sin2^(x2 - ro)), 
(4.23) Vo2 =(0,0,sin2^(a:2 -ro),0), 

^^02 =(cos27r(a;2 - ro), 0,0,0) 
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contribute to the third order terms in evaluation of (|4.21|) . Direct calculation shows 
that 

9fce 



(4.24) 



<G(VL,V'l) 

<G(V'L,V'L) 
<G(vL>V'L) 

<G(V'L,V'L) 



> = 0, 



^-02 > 



^-02 > 



V'02 > 



^l-c 2 



We notice that 

/302 

(4.25) 

/3. 



^-02 > 



<G{4>1,4'1),^02 >=0, 
< 0(^1,^-1), V'^2>= - 

<G(V'L,V'L),V'o2 >=o, 

<G(VL,VL),'^02>=- 
<G«,VL),V'o'2 >=- 
<G(VL,V'feJ,'^02 >=0. 



2 

— ^ ro, 



^fcc 2 



Le ao2 = -47r Le, 



/3o'2 



02 



-Pr"o2 



2Pr 



-47r-^Pr 



rPr. 



By the approximation formula of center manifold functions, see [9], we get 



^/3i,(a;,y) 



(4.26) 



9fc 



SttLc 



Tk 

„2 I „,2 



Plugging ffTFIh into dMU, by and (g^Hl), we obtain 

< G(V'),VL* >= - X < G(V;^^,Vi*),$ > +o(|a;|3 + \y\^) 



(4.27) 



Le 



+ rfe,T^Jx(x2+z/)+o(|a;|3 + |y|3), 



< G(V'), VL* - y < G(V'L,^L*)>* > +oiM' + \y\') 



Te"^" Le 



+ T,Tl)y{x' + y') + o{\x\^ + \y\'). 



By (|4.10p . (|4.13p and (|4.27p . we evaluate (|4?7l) at i? = i?c and obtain the reduction 
equation 



(4.28) 



^^-lxix'+y') + oi\x\' + \yn 



dy 
dt 



-yix' + y') + o{\x\' + \yn 



where b is as defined as 
(4.29) 



6 = 



(l+Lc )aiaTh 



Le3 



(1 + 4;!) + Le^a^ R + Lc f Ii±i^2l£i£Z^ + 
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and aj^ and Rc are given by Remark 13.21 It is obvious that 6 > 0. Hence, standard 
energy estimate gives Assertions (2). Assertions (1) and (3) follows from Theorem 
3.1. This completes the proof of Theorem 4.1. 



5. Convection Scales 

Under the same setting as (|2.17p - (j2.19|) . including the fluid frictions, we consider 
the following nondimensional equation: 

dui ( 2 du2 1 dp\ 
=Prl Aui H — ttUi — OoUi ' 



ot \ To OXi Tq OXi I 



1 

LOU2 ~ (U ■ V)W1 UlU2, 

ro 



— ==Pr( Au2 - — — 2 ""2 + + Rq- diU2 i 



dt \ ro dxi dx2 ' 

(5-1) , , , 1 2 



wui — (u ■ V)u-) H Ui , 

--AT + JU2 + aq-{u- V)T, 



where 



dT 

'dt 

dq 

— =Lc Aq + U2-(u- V)g, 
div u — Q, 



a 

7 = IH , 

' 24' 



(5.2) 



5, = C.,h\ 1 = 0,1, 
Co = 3.78 X 10"^ TO"2 • 
Ci = 6.7 X 10* . s-i. 

The domain is M = [0, 27rro] x (ro, ro + 1), and the boundary conditions are given 

by 

{u,T,q){xi + 27rro,X2) = {u,T,q){xi,X2), 

(5.3) dui 

U2 = 0, - — = 0, T = 0, g = 0, at X2 = ro, ro + 1. 

OX2 

The corresponding eigenvalue problem reads 

^ / . 2 du2 If 9p \ 

Pr Aui H ^ ^0"! - t; — - = pui, 

V ro Oxi rn Ox^ / 



ro r2 dxi ' 

Vv(Au2 - 

(5.4) 



Prf Au2 - A^!iL _ + RT + - 5iU2 - + '^^i = /^"z, 

V ro OTi ri OX2 ' 



AT + 7U2 + = /3T, 
Le Ag + M2 = /3q, 
9mi ^ 9u2 Q 
i9a;2 ' 
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Using the same analysis as in Section 3.1, we obtain the foUowing formula of the 
critical thermal Rayleigh number 

'2 



(5.5) 



where 



i?c = min 

kj>i 



1 

Le 



,4 '0 



kj J^2 ^^k] ^ ^2 



(1 



1 



24 Le aqol^HT^ 



- mm 

k>l 



Le a'i^ X aTarh 
1 + ^. ) 



24 Le aqaf.-^KT 



4 ' 



afeiTi(afei + -2) + -:T + '^i"fei 



*fcl 







fc2 J 


Sqt: 


J 



Let y 



and 
1 



1 



24 



Le aqKT{y + tt^) 



i?+i(y + 7r2)2(y + 7r2 + l) 



+ — (y + TT^) + 5i(?/ + 7r^) + 

^0 

Taking the derivative of g{y), we get 



y 



rg y^ Le a^KTiy + tt^) 



By dO]), we have 

(5.6) Si = 2.76 X 10^°, So = L55 x 10^^ 
Hence, 

(5.7) (5i > (5o > 1. 

Under this condition, the critical value of g{y) is approximated by 

Therefore, the critical thermal Rayleigh number is approximated by 
(5.9) i?c ~ ,5i7r2 = 2.76 X 10^1. 

Next, as in [5J[TT], we adopt 

iy=1.6x IQ-^ m'^/s, kt = 2.25 x 10"^ m'^/s, 
ar = 3.3 x lO^V"*^, Pr = 0.7L 
We notice also that the height of the troposphere is 

/i = 8 X 10^ m. 

Hence, we get 



(5.10) 



(5.11) 



To - Ti = AT, 



Rr = 60° C. 
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Here we notice that the above approximations agree with that of the model of 
tropical atmospheric circulations without humidity. However, as we can see from 
(|5.5p . the coefficient of R is negative. This implies that the humidity factor lowers 
down the critical thermal Rayleigh number a little bit. 

Next, since the nodimensional radius of earth is ro = 6400000//i = 800, we derive 
from 



(5.12) 




that the wave number kc and the convection length scale Lc as 

fcc~6, L = (6400 X 7r)/6 = 3350 fcm. 
This is consistent with the large scale atmospheric circulation over the tropics. 
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